############# added by Steiger
winsorize<-function(x,tr=.2){
#
#  Winsorize a sample
#
#  tr is the amount of Winsorization
#
y<-sort(x)
n<-length(x)
ibot<-floor(tr*n)+1
itop<-length(x)-ibot+1
xbot<-y[ibot]
xtop<-y[itop]
y<-ifelse(y<=xbot,xbot,y)
y<-ifelse(y>=xtop,xtop,y)
return(y)
}

win.mean <- function(x,tr=.2){
win(x,tr)
}

win.var <- function(x, tr=.2){
winvar(x,tr)
}